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4b /n/f/o electronic structure calculations of solid, solution-processed 
metallotetrabenzoporphyrins 

Patrick B. Shea 1 El and Jerzy Kanicki 1 

Organic and Molecular Electronics Laboratory, Department of Electrical Engineering and Computer Science, 
The University of Michigan, Ann Arbor, MI 48109 USA 

An ab initio study of the electronic structures of solid metallotetrabenzoporphyrins (MTBPs) utilized in organic tran- 
sistors and photovoltaics is presented. Bandstmctures, densities of states, and orbitals are calculated for H2, Cu, Ni, and 
Zn core substitutions of the unit cell of solid TBP, as deposited via soluble precursors that are thermally annealed to 
produce polycrystalline, semiconducting thin-films. While the unit cells of the studied MTBPs are nearly isomorphous, 
substitution of the core atoms alters the structure of the bands around the energy bandgap and the composition of the 
densities of states. Cu and Ni core substitutions introduce nearly-dispersionless energy bands near the valence and 
conduction band edges, respectively, that form acceptor or deep generation/recombination states. The following article 
appeared in J. Appl. Phys. Ill, 073709 (2012). 

PACS numbers: 71.20.Rv, 72.80.Le, 81.05.Fb 



I. INTRODUCTION 

Porphyrins serve vital functions in a variety of biochem- 
ical applications, including photosynthesis and catalysis, and 
demonstrate versatile electronic and optical properties. 1 Tetra- 
benzoporphyrin (C36H22N4, TBP) has been demonstrated 
as a semiconductor that exhibits thin-film polycrystallinity 
and field-effect behavior via a soluble precursor molecule, 
wherein a soluble precursor is spun cast to create an amor- 
phous, insulating film, then thermally annealed to form an 
insoluble, polycrystalline, semiconducting thin-film.^— The 
synthetic process used to create the precursor molecules can 
be altered to create variations of the TBP molecule that exhibit 
widely different thin-film morphology and electrical perfor- 
mance; H2, Cu, Ni, and Zn core substitutions (TBP, CuTBP, 
NiTBP, and ZnTBP, respectively) in particular have demon- 
strated organic field-effect transistors (OFETs) with mobili- 
ties (///?£) in the range of 0.01-1 cm 2 /V-Sj^r— TBP has also 
recently been demonstrated as a robust donor material in or- 
ganic photovoltaics (OPVs)»i2iii 

The bulk of theoretical examinations on the electronic 
structures of porphyrins have focused on single molecules, 
typically deriving their interpretation from the Gouterman 
four-orbital model^^ii Many of these single molecule calcu- 
lations have examined the effects of metallation or peripheral 
substitution on the electronic structure of the porphyrins^— 
In the solid state, experimental and theoretical studies on 
the electronic structure and conduction pathways in linearly- 
stacked metallo- porphyrins and phthalocyanines showed in- 
tramolecular charge transport is dependent upon molecular 
spacing, and inter-stack charge transport dependent on molec- 
ular composition^— The prevalence of copper phthalocya- 
nine (CuPc) in OFETs has led to it being studied by ab initio 
calculations both in the single-molecule and solid forms^r— 
Noguchi et al. calculated the crystal growth characteristics 
of solid TBP, reporting that TBP crystallized along the b axis 



parallel to the substrate plane. 36 With the inherent advantages 
of precursor-route metallotetrabenzoporphyrins (MTBPs) in 
OFETs and OPVs, the study and comparison of the electronic 
structures of solid MTBPs is vital to understanding film be- 
havior and device operation. 

In this work, we examine the ground state electronic struc- 
tures of the unit cells of solid, solution-deposited MTBPs uti- 
lized in OFETs and OPVs. Utilizing the CASTEP imple- 
mentation of density functional theory (DFT), bandstmctures, 
densities of states, and orbitals are calculated for MTBPs 
(M=H2, Cu, Ni, and Zn). The solid state bandstmctures are 
compared to the single-molecule results in the literature, and 
analyzed in light of their applications in OFETs and OPVs. 



II. COMPUTATIONAL METHOD 

The theory and implementation of DFT in CASTEP has 
been described by Segall et al£L The calculations described 
here were performed using the spin-polarized local density 
approximation (LDA) with norm-conserving pseudopoten- 
tials to describe the exchange correlation potential. Cal- 
culations were also performed with the Perdew-Wang '91 
(PW91) implementation of the generalized gradient approx- 
imation (GGA) with ultrasoft pseudopotentials and found to 
agree with the LDA results. The planewave cutoff energy and 
k-point set mesh (in a Monhorst-Pack grid) values used in the 
presented results were selected by examining calculation con- 



TABLE I: Parameters for CASTEP calculations. 



^'Electronic mail: pshea@umich.edu http://www.eecs.umich.edu/omelab 



Setting 


Value 


Exchange correlation 


LDA 


Pseudopotential 


Norm-conserving 


Plane-wave basis cutoff 


800 eV 


k-point sampling 


0.02 A" 1 


Scissors operator 


0.0 eV 


DOS smearing 


0.1 eV 



2 




FIG. 1 : (Color online) a) The metallotetrabenzoporphyrin 
(MTBP) molecule and b)the MTBP monoclinic unit cell. 



vergence, and selecting parameter values where changes in 
the total energy (in terms of eV/atom), and in the energy band- 
structure, were negligible with increasing cut-off energy or set 
density. Zero net charge was assumed. The parameters used in 
the calculations are summarized in Table J] The bandstructure 
results are presented as calculated at K, ignoring the band 
smearing and bandgap narrowing that would occur at elevated 
temperature; a Gaussian convolution operator of 0. 1 eV was 
applied to the densities of states to approximate such thermal 
effects. In this work, no scissors operator is applied to the 
bandstructure or densities of states, although LDA is known 
to undercalculate the value of the energy bandgap. 38 The re- 
sults presented are thus compared relative to each other so as 
to discern trends. 

The standard high symmetry notation for primitive mono- 
clinic unit cells was used in this work. 39 The crystal structures 
for the solution-processed MTBPs examined here were deter- 
mined by powder and thin-film X-ray diffraction (XRD), 6-9 
and were found to be nearly isomorphous. All unit cells are 
monoclinic with P2;/n symmetry, with the XRD-derived lat- 
tice parameters reported in Table |II] Relevant XRD-derived 
bond lengths and angles are reported in Table [ill] for the 
solution-processed MTBPs examined here. Also included 
in Table [III] are calculated bond angles and lengths for cor- 
responding single-molecule metalloporphyrins i 22 - 24 ' 40 ^ 3 Prior 
to energy calculations a geometry optimization step was per- 
formed on all unit cells during which the lattice parameters 
were fixed and the atomic coordinates allowed to relax. The 
initial and optimized structures were nearly isomorphous, and 
their corresponding energy band calculations returned results 
with no notable difference. The results presented here are for 
those prior to geometry optimization. 

III. RESULTS AND DISCUSSION 

The bandstructures and partial densities of states (PDOS) 
are displayed in Figures [2] and [3] respectively, with the cor- 



responding, calculated band properties summarized in Ta- 
ble [IV] Presentation of the results proceeds as follows. The 
bandstructure and partial densities of states (PDOS) are dis- 
cussed for metal-free TBP. The analogous results are then 
presented and compared for the metal-substituted TBPs. The 
solid-state results are then compared to the literature on sin- 
gle molecule metallo- porphyrins and phthalocyanines. In 
discussion of materials in the solid-state, the nomenclature 
of valence and conduction bands (VB and CB, respectively) 
is used. When discussing single-molecule results or single 
atoms, the nomenclature of orbitals, such as the highest oc- 
cupied molecular orbitals (HOMO) and lowest unoccupied 
molecular orbitals (LUMO) is used. The top valence bands 
and bottom conduction bands are given particular attention 
due to their direct participation in OFET and OPV operation. 



A. Bandstructure 

The E vs. k and PDOS diagrams for TBP are shown in Figs. 
[2a]and[3a] respectively. Bands appear in pairs due to the pres- 
ence of two distinct molecules in the unit cell, giving rise to 
symmetric and antisymmetric molecular wavefunction combi- 
nations. Degeneracy is noted across much of k-space, except 
between Z-T-Y and B-D-E. The energy bandgap minimum of 
Ec=1.22 eV occurs at k=F. The top valence band, or valence 
band maximum (VBM), is composed of 2p states from the 
porphyrin core C atoms, with a smaller contribution from the 
benzene C atoms. The bottom conduction band, or conduction 
band minimum (CBM), is composed of 2p states from pyrrole 
and meso C and N atoms, and spreads to the C atoms in the 
peripheral benzenes, and are parallel to the molecular plane. 
The top valence band width (W,,) are wider than the bottom 
conduction band widths (W c ) (W v =0.47 eV vs. W c =0. 16 eV), 
although for TBP the bottom conduction bands near the band- 
edge have notably larger densities of states than the bottom 
valence bands. 

The E vs. k and PDOS diagrams for CuTBP are shown 
in Figs. [2b] and [3b] respectively. The VBM and CBM are 
similar to that of TBP, but of interest is the introduction of a 
nearly-dispersionless, degenerate band pair above the VBM 
resulting from the unpaired 3d 9 electron of the Cu 2+ atom in 
a 3d x 2_ v 2 orbital (with significant contribution from the pyr- 
role N atoms), which in single-molecule and solid CuPc is 
reported as the singly-occupied molecular orbital (SOMO). 33 
Thus, here the new band is not designated as the VBM due 



TABLE II: Lattice parameters of reported solid MTBPs 
derived from X-ray diffraction measurements. All cells are 
monoclinic with P2i/n symmetry^— 





TBP 


CuTBP 


NiTBP 


ZnTBP 


a (nm) 


1.2405 


1.2390 


1.2360 


1.2322 


b (nm) 


0.6591 


0.6585 


0.6578 


0.6555 


c (nm) 


1.4927 


1.5150 


1.5190 


1.4994 


ySf) 


101.445 


101.160 


100.620 


101.748 
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TABLE III: Relevant bond characteristics for the examined solid MTBPs derived from X-ray diffraction measurements. Also 
included are calculated bond characteristics of the single -molecule forms of the respective MTBPs i 22,24,42 ' 43 





TBP 


SM TBP 


CuTBP 


SM CuTBP" 


NiTBP 


SM NiTBP b 


ZnTBP 


SM ZnTBP 


M - N (A) 






2.040 


2.029 


1.990 


1.972 


2.089 


2.078 


C9-C10(A) 


1.400 


1.397 


1.383 


1.391 


1.375 


1.381 


1.400 


1.391 


Cll -C16(A) 


1.424 


1.423 


1.415 


1.365 


1.409 


1.362 


1.425 


1.411 


a (M-N-C10) (°) 


123.5 


124.7 


126.4 


124.5 


127.2 


127.8 


124.2 


125.8 


a(C9-C10-Cll)(°) 


126.0 


125.7 


124.0 


126.9 


122.9 


123.7 


125.9 


125.1 


a (N-C10-C9) (°) 


126.9 


126.1 


126.3 


128.1 


126.4 


127.7 


126.9 


125.4 



a Single-molecule form is copper tetraphenlyporphyrin (CuTPP). 
b Single-molecule form is nickel porphyrin (NiP). 





FIG. 2: The energy bandstructure of|(a)]TBP,|(b)]CuTBP,|(c)]NiTBP, and [(d)] ZnTBP. The eV points represent the top of the 

respective valence bands. 



to its composition. The VBM composition is the same as 
TBP, and retains the shape of the valence band as observed 
in TBP. The CBM demonstrates the same orbital composition 
as TBP. The energy bandgap minimum of 1.29 eV occurs at 
k=I\ However, at k=C the states in the CBM are only 20 meV 
higher than for k=r. The difference of less than kT at 300K in- 
dicates that as films are utilized at room temperature, thermal 
effects will narrow the bandgap, smear bands together (see the 
PDOS in Fig. [3bb , and possibly alter the bandgap from direct 
to indirect. The dispersionless SOMO band is a degenerate 
pair with a maximum occurring at k=E, 0.11 eV above the 
valence band; the SOMO to conduction band gap is 1.18 eV. 
At 300K, this degenerate pair is more than kT away from the 
VBM. Like TBP, CuTBP is degenerate except between Z-F- 
Y and B-D-E. The band splitting, however, is less than for 
TBP, particularly in the bottom conduction band. The top va- 



lence bands are wider than the bottom conduction bands, even 
though the bottom two conduction bands overlap (W,,=0.44 
eV vs. W c =0.26 eV). Furthermore, the bottom conduction 
band densities of states have a larger sum of states than the 
top valence bands. 

The E vs. k and PDOS diagrams for NiTBP are shown in 
Figs. [2c] and [3c] respectively. Introduction of the Ni 2+ atom 
with a 3d & orbital induces a nearly-dispersionless, degenerate, 
unoccupied band pair of 3of x 2_ v 2 character, but in this case be- 
low the CBM (also with significant contribution from the pyr- 
role N atoms). As with CuTBP, the VBM and CBM are com- 
posed similarly to TBP. An energy bandgap minimum of 1.35 
eV occurs at k=F. The CBM is shifted upwards relative to the 
valence band, as compared to TBP. The valence band width is 
more narrow than in the other MTBPs (W,,=0.27 eV vs. 0.47 
eV), while the conduction band width is similar to CuTBP 
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(W c =0.25 eV). However, the minimum of the Ni dispersion- 
less band occurs at k=D, with an energy gap to the VBM of 
1.05 eV. Whereas in CuTBP the 3d 9 orbital contributes no- 
tably on the valence band side of the bandgap, in NiTBP the 
Ni 2+ 3d & orbital contributes on both sides of the band gap. 
That is, both the Ni-related mid-gap band, and the band be- 
low the VBM comprise contributions by the Ni 3d & orbital. 
Furthermore, while W v , is more narrow in NiTBP, the PDOS 
reveals that the sum densities of states are notably larger in 
the VBM than for the other MTBPs. As observed in CuTBP, 
degeneracy is more severe for NiTBP than for TBP, with both 
conduction and valence bands displaying smaller band pair 
splitting between Z-T-Y and B-D-E. 

The E vs. k and PDOS diagrams for ZnTBP are shown in 
Figs. [2d]and[3d] respectively. Similar structure is noted in the 
E vs. k diagram compared to TBP, especially in the VBM. 
Band splitting is again observed between Z-r-Y and B-D-E. 
For ZnTBP, however, the band-pair degeneracy in the CBM is 
more pronounced than for TBP between Z-T-Y; the bands are 
also shifted upwards for k^r. The energy bandgap minimum 
of 1.23 occurs at k=F. The VBM is composed is 2p states 
from the porphin core C atoms, and the CBM is composed 
of 2p states from the porphin core C and N atoms. Further- 
more, the contribution of the Zn 3ii 10 orbital near the bandgap 
is confined to a small number of states 1.12 eV below the top 
valence band edge. The valence band width for ZnTBP is 
still wider than the conduction band width (W v ,=0.45 eV vs. 
W c =0.38 eV), although with a notably smaller difference than 
in the other MTBPs. The conduction band width is signif- 
icantly larger than the other MTBPs, although there is little 
notable difference in the magnitude of the PDOS. 



B. Relation to Single Molecule Calculations 

The wealth of literature on single molecule porphyrins and 
phthalocyanines illuminates the bandstructure and PDOS re- 
sults described here for the solid-state. In TabldlVlthe HOMO 
and LUMO orbital energies (Ehomo an d ^lumo-i respec- 
tively) of calculated single-molecule metalloporphyrins are 
presented for comparison to the values calculated for the ex- 
amined solid MTBPs. While the results discussed here im- 
pact charge transport in thin-film OFETs and OPVs, it is dif- 
ficult to speculate on OFET and OPV results at present due to 
the complex nature of those thin films (morphology or contact 
resistance, for example) £~^r— As presented in Table UTR the 
measured bond lengths and angles for the solid MTBPs exam- 
ined here concur with those reported in the single-molecule 
literature, including the trend in M-N bond length. 

The single molecule orbital structure of free-base por- 
phyrins is most commonly associated with the Gouterman 
four-orbital modelr^r— with more recent extensions exam- 
ining various substitutions or phthalocy anines i 3 l i 42A4 ~ £L De- 
termination of the HOMO and LUMO orbitals yields similar 
structure to that observed here in the solid-state. The D21, sym- 
metry of TBP, versus D^, symmetry in the metal-substituted 
molecules, is known to lift the molecular orbital degener- 
acy, as also shown here in the solid-state for k between Z- 



r-Y^s Furthermore, Orti et al. examined the bandstructure of 
one-dimensional stacks of similar H2PC molecules and found 
that notable n-n interactions were exhibited between adjacent 
molecules, but not between molecular columns, similar to the 
solid-state structure here demonstrating weak intermolecular 
interaction^ 

The popularity of CuPc as an organic semiconductor has 
produced many experimental and theoretical studies on the 
electronic structure of its single molecule and solid-state 
formSj 22 i 32 i 33 i 35 i 48-50 T h e aforementioned SOMO level was 
noted in solid CuPc, and the HOMO orbitals similarly as- 
signed to the pyrrole C atoms. 33 The 3d 9 open-shell or- 
bital leaves an unpaired electron in a 3d x i_ Y i orbital with bi g 
symmetry, 14,51 with contributions from the porphin core N 2p 
orbitals. Yang et al. calculated the bandstructure diagram 
for CuPc; the Eg is similar to that reported here for CuTBP, 
although the degeneracy differs due to the slightly different 
symmetry. 35 The minimal variation from single-molecule to 
crystalline form indicates that, as shown here, the intermolec- 
ular interaction is very weak. 32 

In Ni porphyrins, the Ni 2+ atom, with its spin-paired, dia- 
magnetic 3d & open-shell, is found to also introduce a LUMO 
orbital of 3d x 2_ y 2 character . 22,43 ' 52 Krasnikov et al. designated 
the lowest two unoccupied orbitals as LUMO levels, the low- 
est being the Ni 3d x 2 .,2 and the second being related to the 
porphin C and N 2p orbitals. 

Calculations on Zn porphyrins have indicated a closed-shell 
3d 10 orbital structure as displayed here . 14 ' 22 ' 24 Nguyen and 
Pachter examined several ground state zinc porphyrins, in- 
cluding ZnTBP, and reported similar HOMO and LUMO or- 
bital symmetry (2ai H and l eg , respectively)^ The Zn 3d or- 
bital appears in a HOMO-1 orbital (4a2„) approximately 0.7 
eV below the HOMO level; in our bandstructure the Zn 3d 
orbital appears at 1.1 eV below the valence band. Further- 
more, the solid-state bandstructure displays an intermediate 
band composed of 2p orbitals at 0.7 eV below the valence 
band. 

The calculated single-molecule HOMO and LUMO or- 
bital energies from literature shown in Table [IV] illuminate 
many similar trends to those observed here in solid-state 
MTBP calculations. In each case the valence band energy 
(E v ) shifts down in energy relative to the single-molecule 
HOMO energy level. Similar to the calculated properties for 
single-molecule metalloporphyrins, metallated MTBPs dis- 
play a n-n* gap wider than the free-base form. Lower- 
ing of the E v as displayed in the solid MTBPs was also re- 
ported for Ehomo m a series of single-molecule metallated 
tetraphenylporphyrins^ 2 ^ 3 - In both the solid-state and single- 
molecule cases Ni porphyrins display the lowest E v /Ehomo 
and widest band gap. With respect to the mid-gap states in 
CuTBP and NiTBP, the E,,-Es omo separation in CuTBP is 
smaller for solid CuTBP than for the single-molecule forms, 
which suggests a higher degree of Cu-N interaction, while the 
mid gap band-to-E f separation in solid NiTBP is similar to 
that reported in single-molecule nickel porphyrins . 22,43 ' 52 ' 53 
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FIG. 3: (Color online) The densities of states of|(a)]TBP,[(b)]CuTBP,[(c)]NiTBP, and[(d)]ZnTBP. The eV points represent the 

top of the respective valence bands. 



C. Discussion 

The similarities of the energy band structures observed in 
Fig. |2]correspond to the nearly isomorphous crystal structures 
of the examined unit cells. In the MTBP bandstructures re- 
ported here the degeneracy is indicative of the high degree of 
symmetry in the molecules (D^ f° r TBP and for Cu-, Ni- 
, and ZnTBP) and in the monoclinic unit cellM^ It is noted 
in the E vs. k diagrams that the dispersion is very low and 
the bandgap large between k=Y-A and E-C, indicating planes 
of least molecular orbital overlap. In the metal-substituted 
MTBPs examined here the variation in band-splitting may 
be indicative of slightly changed molecularly symmetry or 
planarity, as expected from calculations on single-molecule 
forms of metalloporphyrinsjS21 That the band widths are 
much larger than kT at 300K indicates that only states near the 
VBM and the CBM are likely to be populated. However, in 
less dispersive regions the bands will populate rapidly with in- 
creasing temperature, indicating that the band populations will 
saturate quickly. The trends in W,, and W e imply as well that 
charge will populate the conduction band states, especially in 
less dispersive regions, more rapidly than in the valence band 
states. 

In small-molecule films such as pentacene, transistor mea- 
surements have shown that the herringbone molecular ar- 
rangement, as opposed to a face-to-face arrangement, lim- 
its charge transport by increasing intermolecular spacing and 
reducing n orbital overlap^— Minari et al. examined the 



effects of intermolecular distance and n orbital overlap in a 
series of metal-substituted octaethylporphyrin OFETs with 
the molecules arranged face-to-face along the c axis, and 
found fif£ to be inversely dependent on the intermolecu- 
lar spacing. 58 XRD measurements on MTBP thin-films have 
shown that molecules stack in the herringbone arrangement 
along the b axis, with the molecule face parallel to (1 14) and 
the HOMO orbital surface parallel to the molecule face. 36 
Thus, in a and c the intermolecular spacing is small, but n 
orbital overlap is also small. In b, however, the intermolecular 
spacing is small, and the n orbital overlap is enhanced. 

That b is the preferred path of charge transport is confirmed 
by examining the curvature of the band diagrams in Fig. |2] 
The curvature of the VBM between F-Y indicates the direc- 
tion of the smallest hole effective mass (m?). Conversely, 
the valence band splitting between T-Z (along c) indicates the 
lower band displays a positive ml, a reflection of the herring- 
bone alignment of the molecules. It is also worth noting that 
the curvature of the VBM in NiTBP corresponds to the lesser 
W,,, such that m* h for NiTBP would be expected to be larger 
than for the other MTBPs. Furthermore, the curvature of the 
conduction bands is lesser than for the valence bands near the 
band gap, indicating that the electron effective mass (m*) is 
larger than m^. 

While the exact nature of the effect of the Cu and Ni 
3flf x 2_ v 2 orbitals on solid, thin-film charge transport is out- 
side the scope of the analysis performed here, for example 
due to the widely varying aggregate nature demonstrated in 
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TABLE IV: Calculated energy bandgaps (Eg), valence band energy (E,,, relative to the vacuum level), valence band widths 
(W v ), and conduction band widths (W c ) of solid MTBPs. Also included are relevant, calculated single-molecule (SM) HOMO 

and LUMO energy levels . 24 ' 42 ! 53 



Material 


Eg (eV) 


Ev (eV) 


W,. (eV) 


W c (eV) 


SM Ehomo (eV) 


SM Elumo (eV) 


TBP 


1.22 


-5.11 


0.47 


0.16 


-4.77 


-2.32 


CuTBP^ 


1.29 


-4.96 


0.44 


0.26 


-4.89 


-2.37 


NiTBP^ 


1.35 


-5.49 


0.27 


0.25 


-4.96 


-2.30 


ZnTBP 


1.23 


-5.42 


0.46 


0.38 


-4.63 


-2.11 



a Single-molecule form is CuTPP. 
b Single-molecule form is NiTPP. 



solution-processed MTBP films, 4 ' 7-9 given the relative results 
presented here the nature of those bands could be expected to 
affect charge transport as follows. The curvature of the dis- 
persionless, mid-gap bands induced by the Cu 2+ 3d 9 and Ni 2+ 
3c/ 8 orbitals in the 3d x 2- y 2 configuration indicates mjj in the 
CuTBP band is large, and an electron excited to the NiTBP 
band will be expected to have a large m* compared to the 
CBM. The partially-filled Cu bands, with rising temperature, 
accept an electron from the valence band and hence act as a 
deep acceptor; from the PDOS diagram this acceptor would be 
expected to contribute charge carriers in an OFET or OPV. In 
NiTBP, there are two bands of interest. The unoccupied mid- 
gap band, given that it is 0.3 eV from the CBM, potentially 
behaves as a deep electron trap. The band below the VBM, on 
the other hand, would be expected to contribute to charge con- 
duction in an OFET or OPV. This could potentially explain 
the higher fj.fg noted in OFETs fabricated from CuTBP and 
NiTBP (>0.1 cm 2 /V-s) compared to TBP and ZnTBP (0.01 
cm 2 /V-s), although the presented calculations cannot account 
for grain boundary regions observed in the widely-varying 
thin-film aggregation/ 4 * 7 ^— 



Charge transport measurements and bandstmcture calcula- 
tions on one-dimensional stacks of metallo- porphyrins and 
phthalocyanines, including H2, Ni, and Zn complexes, have 
indicated that the interaction of the n orbitals from the por- 
phyrin ligand with the metal d x 2_ y 2 orbital can play a signifi- 
cant role in intramolecular orbital formation and intermolec- 
ular charge transport along the one-dimensional stack and in 
bulk crystals^— Ishikawa, for example, found by calcula- 
tion and solution absorption spectroscopy that linear stacks 
of closed-shell phthalocyanines could be treated as metal- 
free phthalocyanines. 29 Comparatively, whereas conduction in 
one-dimensional stacks of NiPc has been suggested to result 
solely from ligand orbitals, 59 conduction in one-dimensional 
NiTBP has been shown to involve both the ligand n orbitals 
and the metal d orbitals. 60 The limiting factor in charge trans- 
port within the unit cell can thus be expected to be the n or- 
bital overlap resulting from the herringbone molecular pack- 
ing. Limited enhancement to charge transport can be ex- 
pected from 3d orbital introduction, even with the introduc- 
tion of charge-contributing acceptors, without improvement to 
molecular packing, such as in the [i^e enhancement in OFETs 
observed in TIPS-pentacene compared to pentacene i 57 ' 58 ' 61 



IV. CONCLUSIONS 

The energy bandstructures and densities of states of 
solid metallotetrabenzoporphyrins used in solution-processed 
OFETs and OPVs have been investigated and analyzed. Cal- 
culations indicate the effect of metallation depends on the or- 
bital structure of the substituent, such that metals with open- 
shell 3d orbitals (Cu and Ni) induce interactions between 
the 3d states of the metal and n states of the tetrabenzopor- 
phyrin molecule near the bandgap, potentially acting as ac- 
ceptor states in Cu and deep trap states in Ni. An empty core 
(H2) or one with a full 3d orbital (Zn) are similar, with neg- 
ligible contribution of d states near the bandgap. Otherwise 
the structures of the unit cells are nearly isomorphous and i so- 
electronic. Furthermore, the bandstmcture diagrams indicate 
enhancement is limited by weak intermolecular interaction. 
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